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Flat directions in the minimal supersymmetric standard model are known to deform into non- 
topological solitons, Q-balls, which generally possess both baryon and lepton asymmetries. We 
investigate how Q-balls evolve if some of the constituent fields of the flat direction decay into light 
species. It is found that the Q-balls takes a new configuration whose energy per charge slightly 
increases due to the decay. Specifically, we show that all the stable Q-balls eventually transform 
into pure B-balls via the decay into neutrinos. 
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I. INTRODUCTION 



Non-topological solitons, such as Q-balls 1|] and I- 
balls 01, play important roles in the particle cosmology, 
partly because they are (quasi-)stable objects. Their sta- 
bility comes from the conservation law. The Q-ball is 
composed of a complex scalar field with an U{1) charge, 
the conservation of which enables the Q-ball to be sta- 
ble or long-lived. Similarly, /-balls are composed of a 
real scalar field (or a complex scalar field in a nearly 
straight-line motion), whose dynamics conserve the adi- 
abatic charge. Since we can only see either these objects 
in their final state or their decay products in the present 
universe, it is important to study their evolution during 
the course of the expanding universe. 

The baryon asymmetry is one of the key observational 
results to uncover the history of the universe. Among 
many baryogenesis scenarios proposed so far, the mech- 
anism proposed by Affleck and Dine 3] can be realized 
in the minimal extension of the standard model, that is, 
the minimal supersymmetric standard model (MSSM). 
The mechanism makes use of one of flat directions, along 
which there is no classical potential in the supersymmet- 
ric limit. Although the fiat direction is parametrized by 
a set of chiral superfields, it is convenient and even suf- 
ficient in the usual situation to express it in terms of 
a single complex scalar field dubbed 'Affleck-Dine (AD) 
field,' <l>, which generally has nonzero baryon and lepton 
charges. The dynamics of the AD field not only gen- 
erates the desired baryon asymmetry, but allows non- 
topological solitons to be formed. The properties of Q- 
balls in the context of the AD mechanism have been ex- 
tensively studied 0, 0, IE 0- Q-balls can be stable or 
unstable depending on the situation. Stable Q-balls can 
be dark matter 0^0, 0, which might solve the coinci- 
dence problem between the abundance of baryons and 
that of dark matter. Unstable but long-lived Q-balls can 
keep the baryon and lepton asymmetries stored inside 
them from being subject to the sphaleron effects. For in- 
stance, they can successfully generate both large lepton 
asymmetry and small baryon asymmetry 10]. It is even 
possible to realize the late-time baryon appearance after 
the relevant BBN epoch with use of Q-balls |Tl|. /-balls 
are also formed in the AD mechanism as an excited state 
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of Q-balls. 

In the examples mentioned above, the single-field 
parametrization is simple and valid. However, if we ex- 
amine the evolution of Q-balls closely, we have to pay 
much more attention to the fact that the flat direction 
is actually composed of several scalar fields. Indeed, the 
slepton condensate always decays into a pair of neutri- 
nos, while the squark condensate can decay into hadrons 
only if its energy per one baryon charge exceeds the nu- 
cleon or meson masses. Therefore, if a Q-ball has nonzero 
baryon and lepton asymmetries and its energy per unit 
charge is small enough, the leptonic component decays, 
leaving only the baryon asymmetry inside Q-balls (i.e., 
pure 5-balls). However, since the Q-ball solution is ob- 
tained under the assumption that the U{1) charge, Q, 
is conserved, it is not trivial how the Q-ball configura- 
tion changes during the course of the partial decay. Do 
Q-balls come apart? Or do they change their shape to 
satisfy the modified Q-ball solution? It is the purpose of 
this study to answer these questions. 

In this paper, we investigate how Q-balls evolve if some 
of the constituent fields decay into lighter species. In the 
next section, we review the flat direction of the MSSM 
and the Q-ball solution. In section IhTI we study the par- 
tial decay of Q-balls with a toy model of the flat direction. 
The results obtained there are confirmed with use of nu- 
merical calculation presented in section IIVI Finally we 
give discussions and conclusions in section IVl 



II. PRELIMINARIES 

In this section, first we explain how fiat directions are 
expressed, paying particular attention to the single-field 
parametrization. Then we review the Q-ball solution rel- 
evant for the next section. 



A. Flat directions 

In the MSSM, there are many flat directions along 
which both the D-ieiia and F-term potentials vanish at 
the classical level. The /^-flat direction is labeled by a 
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holomorphic gauge- invariant monomial, X. Most of D- 
flat directions can be also F-flat, simply due to a gen- 
eration structure of the quark and lepton sector. Since 
it is easy to satisfy the F-flat condition, let us concen- 
trate on the D-flat condition in the following. The D-flat 
direction, X, can be expressed as 



N 



(1) 



where N superfields {^i} constitute the flat direction 
X, and we have suppressed the gauge and family indices 
with the understanding that Latin letter i contains all 
the information to label those constituents. When X has 
a non-zero expectation value, each constituent field also 
takes a non-zero expectation value 



(2) 



Here each / \/2 is the absolute value of the expectation 
value and is related to every other due to the D-fiat 
conditions: 



0, 



(3) 



where T4 are hermitian matrices representing the gener- 
ators of the gauge algebra, and they are labeled with A. 
This condition can be usually satisfied if we take all 
equal: = Note that D-fiat condition dictates 

that all the amplitudes be equal, while the phases, 6^^, 
remain to be arbitrary. In fact, it is necessary to know 
the pattern of baryon and lepton symmetry breaking, in 
order to identify the Nambu-Goldstone (NG) boson rel- 
evant for the AD mechanism. The situation becomes 
simple if the spontaneously broken symmetry coincides 
with the explicitly violated one. This is the case if the 
A-term is given by some powers of X. Then all Oi become 
equivalent, and the NG boson can be identified with the 
average of Oi [l^] : 







(4) 



where is orgthogonal to the NG mode corresponding 
to the spontaneously broken U{1) gauge symmetry, that 
is, the linear combination of left isospin Tls and weak 
hypercharge Y. This is because TLs,i = = 

by definition of the D-fiat direction. Thus, the dynamics 
of the fiat direction, X, can be well described by one 
complex scalar field, <l>, which is defined as 



V2 



(5) 



Note that ^ is canonically normalized. 

Flat directions are lifted by the supersymmetry 
(SUSY) breaking effect and non-renormalizable terms. 



Since we are interested in the evolution of Q-balls which 
are formed after the fiat direction starts to oscillate, we 
neglect the non-renormalizable terms in the following 
discussion. To be concrete, let us adopt the gravity- 
mediated SUSY breaking model. The potential of the 
fiat direction is now given by 



Vim) = E"^ii**ni+^*iog 



A \ij J 
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(6) 



where rui is a soft mass, Ki a coefficient of the one- 
loop correction, the renormalization scale to define 
the mass. This potential reduces to the following in the 
single-field representation. 



1 + Klog 



where and K are defined as 



1 

N 



E 



K 



(7) 



(8) 



As long as the Z)-fiat condition is satisfied and all 
can be treated equally, the single-field parametrization is 
useful, and we only have to deal with the simple potential 
Eq. {TJ. However, once this assumption breaks down 
(say, some of the constituent fields decay), it is necessary 
to adopt the full potential Eq. Q in order to follow the 
evolution of each . 



B. Q-ball solution 

During infiation, the fiat direction, X, is assumed to 
take a large expectation value (X) 7^ 0. After inflation 
ends, the Hubble parameter starts to decrease, and be- 
comes comparable to the mass of ^ at some point. Then 
^ starts to oscillate in the potential given as Eq. Q. 
However, if the sign of K is negative, it experiences spa- 
tial instabilities, leading to the Q-ball formation 
Here we write down the equation that dictates the con- 
figuration of the Q-ball and derive the solution of the 
gravity-mediation type Q-ball. 

The Q-ball configuration is the solution that minimizes 
the energy E at fixed U{1) charge Q, where E and Q are 
defined as 



E 



I d^x|^(|^,|^ + |V$,|^)+n{^J) 



(9) 
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In the single-field parametrization, they become 
E = j (fx^\^\^ ^\V^^ ^V{^)Y 



(10) 



Here qi and q are the U{1) charges of and respec- 
tively, and they are related to each other as 



N ^ 



(11) 



With the use of the method of Lagrange multipliers, the 
problem is reduced to minimizing 



= y"d3a;{|$|2 + |V$|2 + y($)} 
+w Q + i J (fx q (^$*$ - 4*$ 



, (12) 



where a; is a Lagrange multiplier. The solution takes 
the spherically symmetric form <l>(r, = 0(r) e^'^'^^/v^, 
where ^(r) satisfies the following equation. 



r dr 

with the boundary conditions 



dr 



0, 
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0. 



(13) 



(14) 



r=0 



Here we defined K;(^) = ^(^) — q^oo'^cj)^ /2. For the solu- 
tion to exist, V{(j))/(jP' must have a global minimum away 
from the origin. 

Let us consider the case that the potential is given by 
Eq. Q). The Q-bah equation, Eq. then reads 



ar^ r dr 

-m%(}) ( 1 + K + Klog 



0. (15) 



If we adopt the Gaussian ansatz, 0(r) = 0(O)e~^ ^^Q, 
the radius and angular velocity of the gravity-mediation 
type Q-ball are determined as 



Rc 



V2 



-, qcj = m^\/l + 2\K\. 



(16) 



where we set y]VM* = 0(0)/ V^. The energy per unit 
charge, E/Q^ can be calculated with use of this solution: 



I = ^(l+^w)(l+2|i^|)-/^ 



(17) 



where we assumed \K\ <C 1. Note that these results 
depend on the dimension of space, D = 3. For later use, 
we write down the result in the case of D = 1: 



Rq = 



E 
Q 




(18) 



III. MODEL 

Now let us consider the evolution of Q-balls, assum- 
ing that some of the constituent fields decay into some- 
thing else. For a concrete discussion, we take up the 
simplest possible flat direction composed of two scalar 
fields with U{1) gauge symmetry, X = <l>i(+l) <l>2(— 1), 
but it is trivial to extend our results to more generic 
case. We take the global U{1) charges of <l>i and ^2 are 
Qi = Q2 = Q = ^ ^- Also we assume that only $2 decay 
into light species. After X starts to oscillate, it feels spa- 
tial instabilities and deforms into Q-balls. Once Q-balls 
are formed, the evolution of the scalar fields inside them 
decouple from the cosmic expansion, so we neglect the 
effect of the expanding universe in what follows. 

The constituent fields, <^>i,2, obey the equations of mo- 
tion ^, 



^2 + r$2 + U = 



(19) 



with 



m,*2) = X^mf|$ip(^l + ifaog(K 



(20) 



where F is the decay rate of $2, ^ the U{1) gauge coupling 
constant. Before the decay becomes effective, the Q-ball 
solution takes the form 

$W = $W = * = Me-*-V< (21) 
V 2 2 



^ Even if $1 and ^2 are charged under different global U(l) 
symmetries (e.g., the baryon and lepton symmetries), the sys- 
tem is reduced to that with single global U(l) symmetry (e.g., 
U{1)b-\-l), as long as the A-term is given by some powers of X. 

^ The decay of the Q-ball occurs only around its surface due to the 
Pauli blocking 13j, if the decay particles are fermions. Therefore 
the empirical treatment of inserting r<i>2 might not be valid in the 
light of microphysical processes. However, what we are concerned 
here is not the decay process but the final state after the decay 
process completed. Therefore we believe such a simplification 
does not spoil our discussions. 
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with Rq and uj given in Eq. (|TB|) (or (|T^ ). and 0(0) = 

If we observe the I^-term interaction between 
and $2, the final state can be inferred by assuming that 
the Q-bah remains to be stable during the decay, which 
will turn out to be valid later. First of all, since only the 
absolute values of ^1,2 participate in the D-teim inter- 
action, they can exchange the energy but not the global 
charge. Second, the energy of cannot be extracted 
without extracting its charge, because the motion of <l>i 
is a circular orbit around the origin. Therefore, will 
continue moving in a circular orbit, while the angular 
velocity of <l>2 decreases to zero with the I)-fiat condi- 
tion, 1^2! = 1^1 1, satisfied. The final state can be thus 
expressed as 



$(/) ^ |$(/)| = Me-'-Vi?^ 

2 I 1 I 2 



(22) 



Note that the amplitude, radius, and angular velocity, 
0(0), Rq^ and a), do not necessarily coincide with 0(0), 
i?Q, and cj, respectively. 

Let us derive the Q-ball solution in the final state. 
From the above discussion, we would like to find the so- 
lution of the following form. 



$1 = — , ^2 = |^l| 



(23) 



with i>2 = 0. It should be noted that ^2 still has the 
gradient and potential energies in order to satisfy the D- 
fiat condition. The energy and charge of the system are 
given as 



E 



$1 



lv*l' 



+m||$|2 {l + K\oi 



Q = 



where we defined 



mf + mi 



K 



m{Ki + m2i^2 
m\ + m2 



As usual, we find the solution that minimizes 



Q + ^ j (fx (<i>*i> - <i>*<i>) 



(24) 



(25) 



(26) 



where w is the Lagrange multiplier. The solution is given 
as $(r, = (^(r)e'"*/v^' where 0(r) satisfies 



cP4> ^ 2d(j) ^ u!^^ 
dr^ r dr 2 

-m%cj)(l + K + K\og( 



V4M2 



(27) 




0, 



FIG. 1: The schema of the effective potential of $i . Before $2 
decays, $1 inside the Q-ball moves in a circular orbit shown 
as the dotted line. After the decay of $2, the amplitude of $1 
decreases but its angular velocity increases so that the charge 
of $1 remains constant. The modified orbit is represented as 
the solid circle. 



with the boundary condition given as Eq. (|14|) . Assuming 
the Gaussian ansatz, 0(r) = 0(O)e~^ /^q^ we obtain 




log 



^(0) 
0(0) 



, (28) 



Meanwhile, the charge conservation of <I>i leads to 

cj0(O)^ =a;0(O)^ (29) 

Assuming \K\ <C 1, co and 0(0) can be iteratively deter- 
mined as 



CO 

m 



(2 + 4|if| + |i^|log2)i/2m$, 

(2 + 4|if| + |if|log2)-i/V(0), (30) 



The energy per unit charge, E/Q, can be calculated with 
use of this solution, and given by 



E 



V2m$ (^l + ^|if| + ^log2y (31) 



Thus the Q-ball configuration after the decay of <I>2 has 
the same shape in the real space, but the scalar field 
inside the Q-ball moves in a different orbit. We write 
down the similar results in the case of D = 1. 



Rq 



V2 



(2+|if|log2)^/^m$, 

1/4 



m - {2+\K\\og2)-''''m, 



(32) 



We would like to make several remarks. First, the Q- 
ball solution such as Eq. exists if and only if the 
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Q-ball solution is allowed before the decay of ^2, that 
is to say, K < 0. In other words, the partial decay of 
the constituent fields does not split the Q-ball into frag- 
ments. This is rather surprising result. For instance, 
even if Ki is positive, the Q-ball configuration remains 
stable as long as K is negative. It is the conservation of 
the charge of <I>i that keeps the Q-ball stable during the 
course of the decay. After the decay of ^2 completed, $1 
still feels the potential of ^2 through the D-teim poten- 
tial. Therefore moves in an effective potential given 
as the sum of the potentials of and ^2- See Fig. ^ 
Second, the energy per unit charge becomes about >/2 
times larger, which means the Q-ball configuration be- 
comes slightly unstable. The resultant Q-ball might be 
able to decay into light particles by using this increase. 
Such increment comes from the fact that the interplay 
of the charge conservation of ^1 and the D-fiat condi- 
tion do not allow ^2 to decay with the potential energy. 
Third, the reason why the energy of cannot dissipate 
through the decay channel of $2 is that the Q-ball is a 
state of the lowest energy. It is helpful to consider how 
the /-ball 0, which is considered to be in an excited 
state of the Q-ball, change its profile if the partial decay 
occurs. The /-ball solution exists if the scalar potential is 
dominated by the quadratic term, which guarantees the 
approximate invariance of the adiabatic charge, /. The 
/-ball configuration is obtained as a state with the low- 
est energy at fixed /. Its profile is very similar to that of 
the Q-ball, except that the scalar field moves in an orbit 
with large ellipticity. If the kick of A-term is not strong 
enough, /-balls are generally formed instead of Q-balls. 
Let us consider the /-ball composed of X = ^1^2, and 
suppose that ^2 decays into something else. In this case, 
the energy of <l>i, in addition to $2, can dissipate into 
light particles through the D-teim interaction. This pro- 
cess lasts until the ellipticity of the orbit becomes ~ 1, 
i.e., until the /-ball becomes a Q-ball, from which the 
energy cannot be extracted anymore without violating 
the charge conservation of ^1. 



IV. NUMERICAL CALCULATION 

We have performed the numerical calculation to vali- 
date the analysis presented in the previous section. In 
particular, it is important to follow the evolution and 
show that the Q-ball does not break up during the decay 

of <I>2. 

We have solved the equations of motion, Eq. (|T^ . on 
the one dimensional lattices with the initial condition 
given as Eq. (|21|) . The space x and time t are normalized 
by a mass scale m, while the field value is normalized by 
M*. We take the following values for the model param- 
eters: mi = m2 = m, (/)(0)/M* = 2, g'^M'l = 200m^, 
Ki = 0.02, K2 = -0.1, and F = 0.1m. 

The evolutions of |^i| and |<I>2| are illustrated in Fig. [21 
As seen in the figure, the /^-fiat condition is satisfied in- 
side the Q-ball but not outside. However, as shown in 
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FIG. 2: Theevolutionofthedistributionof |$i| and |$2|. The 
horizontal straight line represents the analytical estimate of 
(j) given in Eq. (|32|) . The /)-flat condition is satisfied inside 
the Q-ball centered at x = 100. The span of the lattices is 
actually x = - 200. 

Figs. El and 21 the leakage of the energy or charge of <I>i 
from the Q-ball is negligible ^. The amplitude decreases 
as expected by the analytical estimate (see Eq. ^^). 
The charge and energy of ^i, Qi and Ei, do not change 
before and after the decay. See Figs. 01 and ^ On the 
other hand, the charge of <I>2 quickly decreases to zero, 
and its energy, E2, is also brought down to the potential 
energy of ^2, which is necessary to satisfy the /^-flat con- 
dition. It should be noted that the shape of the Q-ball re- 
mains unchanged throughout the decay as expected. The 
behavior of the energy per unit charge, E/Q, is shown in 
Fig. and it agrees well with the predicted values given 
in Eqs. (tTH|) and 



V. DISCUSSIONS 

So far we have presented the analytical and numeri- 
cal study of the partially decaying Q-balls with use of 
a toy model of fiat directions. Let us consider now the 



^ For a long-time simulation, the leakage of the charge and energy 
from the Q-ball cannot be negligible, as suggested by the slight 
decrease of Qi and Ei in Figs. |3 and ^ respectively. In the real- 
istic situation, the D-flat condition is much severer than in the 
numerical calculation, which would highly suppress such leakage. 
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FIG. 3: The evolution of the charge distribution of $i and 
$2. While Qi is conserved, Q2 rapidly decreases to zero due 
to the decay of $2- The charge outside the Q-ball is negligibly 
small. 
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FIG. 4: The evolution of the energy distribution. While Ei 
stays unchanged, E2 decreases to some finite value, which 
corresponds to the potential energy of $2- We do not have 
included the D-term potential in the definition of ^1,2- 
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FIG. 5: The time evolution of the energy per unit charge, 
E/Q. The analytical estimates for the initial and final states 
are shown as the horizontal straight lines. 



implication of the results presented above. We have as- 
sumed that only a part of constituent fields decay, but 
our arguments apply to more generic case where there is 
a hierarchy between the decay rates, say, Fi ^ F2. The 
asymmetry of the Q-ball continuously changes during the 
course of the decay. Most of the flat directions in the 
MSSM are composed of both squarks and sleptons, there- 
fore they have both baryon and lepton number. If such 
a flat direction deforms into Q-balls and they are stable, 
our argument can be applied since the slepton conden- 
sate can always decay into a pair of neutrinos through 
the exchange of the gauginos. All through the decay 
processes, the Q-balls sequentially transform themselves 
toward pure 5-balls with no lepton number. In the usual 
scenario considered so far, it is a flat direction composed 
only of squarks that forms pure 5-balls. Therefore our 
analysis indicates that the pure 5-balls are generated for 
most of flat directions, as a result of the decay into neu- 
trinos. Specifically, all the Q-ball dark matter must be 
pure 5-balls. Interestingly enough, the resultant B-hal\ 
might become unstable due to the increase of its energy 
per unit charge, leading to a complete decay of the Q-ball. 
In connection with the transformation of Q-balls into B- 
balls, the leptogenesis might become possible even for a 
flat direction with B — L = 0. Another interesting possi- 
bility is that the /-ball can be naturally altered to Q-ball, 
if one of the constituent fields decay into something else. 
In particular, the /-ball composed of a flat direction that 
includes a leptonic field is generally transformed into a 
Q-ball. In fact, since the /-ball is a quasi-stable object, 
its conversion to a Q-ball was not known. The results of 
this study affords some new perspectives on the decay or 
stabilization processes of non-topological solitons made 
of flat directions. 

In summary, we have investigated the evolution of Q- 
balls, assuming that some of the constituent scalar fields 
decay into light particles. With a toy model of the fiat 
direction, we have obtained the analytical solution of the 
Q-ball in the final state, and found that the spatial shape 
of the Q-ball does not change, but the orbit of the scalar 
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field inside the Q-ball does ciiange in order to compen- 
sate for the decay. In particular, the energy per unit 
charge generally increases due to the partial decay, which 
might induce further decay of the remnant scalar fields. 
Also we have performed numerical calculations to con- 
firm these results. Again we stress that Q-balls remain 
Q-balls throughout the partial decay, due to both the 
charge conservation of the remnant scalar fields and the 



i^-fiat condition. 
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